CORRIGENDUM 



Two-dimensional electron gas in a uniform magnetic field in the presence of a <5-impurity 

R M Cavalcanti and C A A de Carvalho 1998 J. Phys. A: Math. Gen. 31 2391-9 



Footnote 1 (p. 4 of the e-print :ond-mat/9704156, p. 2393 of the published version) propagates an erroneus as- 
sertion made in jjj, namely that "in dimensions D > 2 a repulsive ^-interaction expels the S- waves from the Hilbert 
space." By treating a repulsive delta-function potential as the limit of a sequence of square barriers, one can show 
instead that in two or more dimensions the iS-waves are the same as those in the absence of the potential, except at 
the position of the potential, where they drop to zero discontinuously Hence the repulsive 5-potential is trivial 

— it does not scatter and does not affect the energy levels. 

None of the results or conclusions of our paper is affected by this correction. 
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Abstract 



The density of states and the Hall conductivity of a two-dimensional elec- 
tron gas in a uniform magnetic field and in the presence of a (5-impurity are 
exactly calculated using elementary field theoretic techniques. The impurity 
creates one localized state per Landau level, but the Hall conductivity is unaf- 
fected. Our treatment is explicitly gauge- invariant, and can be easily adapted 
to other problems involving zero-range potentials. 
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I. INTRODUCTION 



One of the most puzzling features of the quantum Hall effect is the apparent insensitivity 
of the quantization of the Hall conductivity an (in multiples of e 2 /h) with respect to type 
of host material, geometry of sample, presence of impurities, etc. Prange was probably 
the first to address the question of the influence of impurities on the quantization of <jh- 
For a two-dimensional electron gas in crossed electric and magnetic fields in the presence 
of a 5- impurity, he showed that a localized state exists, which carries no current, while the 
remaining nonlocalized states carry an extra Hall current which exactly compensates for the 
part not carried by the localized state. 

Notwithstanding Prange's claim that his calculation was exact, he had in fact to resort 
to some approximations. In part, this occurred because he worked with a finite (although 
small) electric field, and in part because he used a 5-function potential, which is too 'strong' 
in two (or more) dimensions |JJ, even in the presence of a magnetic field. As shown explicitly 
by Perez and Coutinho Q, if one solves the Schrodinger equation for a square well of radius 
a and depth Vq{o) ~ a~ 2 , one finds a bound state whose energy — oo when a — > 0. In 

order for Eb to remain finite in the limit a — > 0, the depth of the well must diverge more 
slowly than a~ 2 ; explicitly, Vo(a) ~ l/a 2 \n(a/R) (R is some constant with dimension of 
length). Working with this regularized version of the 5- function, they managed to find a 
spectrum similar to the one found by Prange. 

The purpose of this paper is to revisit this problem using elementary field theoretic 
techniques. There are a couple of reasons for doing things this way: (i) it is very easy to 
find the Feynman propagator in the presence of a 5- impurity |5|||, and (ii) the density of 
states and the conductivity tensor can be computed exactly, in an explicitly gauge-invariant 
way. The very singular nature of the 5-function potential in two dimensions shows up in our 
treatment as an infinity in the propagator, but to deal with it is a very simple exercise in 
renormalization ||. 

Our calculations essentially confirm Prange's results. 

II. THE FEYNMAN PROPAGATOR 

Let us consider an electron gas in two dimensions in a uniform magnetic field, in the 
presence of a (^-function potential at the origin. Its lagrangian density is given by (we use 
units such that m = H = c = 1) 

C = ^(idt-H + fi)i). (1) 
The '1-particle hamiltonian' H can be split in two pieces: H = H + V(x), where 

H = - ieA) 2 , V(x) = A5 2 (x). (2) 

The vector potential A generates a uniform magnetic field B (= <9iA 2 — c^Ai), and \i denotes 
the chemical potential. 
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As we shall see in the next section, the particle density and the conductivity of the system 
can be computed once one knows the Feynman propagator, which satisfies the following 
equation (x = (t, x)): 

(id t -H + n) x G(x,x') = 5 3 (x-x'). (3) 
Since H is time-independent, we can look for a solution of @ in the form 

G(x,x')= r^e-^')G(c;x,xO. (4) 

J-oo 2,71 

This in fact is a solution of (H), provided G(u; x, x') is a solution of 

{u -H + n) x G(u;x,x!) = 5 2 (x-x). (5) 
This can be solved in the usual way as 

*xy) = £«^), (6) 

where E n and ^n(x) are the eigenvalues and eigenfunctions of H, respectively Since if is a 
hermitian operator, its eigenvalues are real, so that a prescription must be provided to deal 
with the poles of G(u; x, x') when one performs the integral over u in (f|). For the Feynman 
propagator, this amounts |7|] to deform the integration contour in the complex cu-plane as 
indicated in Figure 1. 

Defining the 'unperturbed' propagator G (uj; x, x') as the solution of (|5|) with V = 0, we 
can formally solve for G as 

G = G + GqVGq + G VG VG + ... (7) 

(The product of two operators F and G is defined as (FG) (x, x') = / d 2 yF(-x, y) G(y, x'), 
and V(x, x') stands for y(x)5 2 (x — x').) Because of the very simple form of V, all the 
integrals in (0) can be performed exactly, and the series can be summed in closed form (we 
drop the dependence on to for simplicity): 

oo 

G(x, x') = G (x, x') + A G (x, 0) £ [A G (0, 0)] n G (0, x') 

Go(x,0)Go(0,xQ 
" G0(X ' X)+ i-G (0,0) ' (8) 

One can verify, by direct substitution in Eq. (]5|) (with V = 0), that the 'unperturbed' 
Feynman propagator is given by (eB is assumed positive) 

G (u; x, x') = — M(x, x') e- B ( x - x ') 2 / 4 V Me£ (x - x') 2 /2) 

2tt V ; ^ u - (n + 1/2) eB + ^ KJ 

where L n (z) is a Laguerre polynomial || and M(x, x') is a gauge dependent factor, which 
can be written in a gauge covariant way as 
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M(x,x') = exp lie J f A(z) • dzj 



(10) 



the integral being performed along a straight line connecting x' to x. 

Given the explicit form of G (u; x, x'), (§) gives the solution of Eq. (||), but, as it stands, it 
is meaningless: the denominator of the second term on the r.h.s. is logarithmically divergent. 
However, this divergence can be absorbed in a redefinition of the 'coupling constant' A: 
introducing a convergence factor e~ an in the sum over Landau levels, one finds (2 = 1/2 — 
(u + n)/eB) 



eB 00 

G (^;0,0) = — £ 



2vr ^ u - (n + 1/2) eB + [i 



1 

2^ 



Ln=0 



z + n n + 1 
(a -> 



E 

n=0 



n + 1 



« — [7 + VK-s) + In a] 

Z7T 

t/>(2) denotes the digamma function and 7 = 0.577 ... is Euler's constant 
a renormalized 'coupling constant' \r as 

1 1 1 / 1 \ 
— = -- — (7 + lna), 

AR A Z7T 



(11) 

Now, we define 
(12) 



and make A depend on a in such a way that A# remains finite^] in the limit a — > + . In 
terms of the renormalized 'coupling constant,' Eq. (P) becomes 



G (x,0) G (0,x') 



G(x,x') = G (x,x') + T _ ^ 



eB 



Go(x,x') + Gi(x, x') 



(13) 



which is now well defined. 



III. PARTICLE DENSITY, LOCALIZED STATES AND HALL CONDUCTIVITY 



A. Particle density 



The particle density is given by 



nix) = — i lim G(t,x.:t' ,x). 

v ' f->t+o 



(14) 



The 'unperturbed' part of n(x) is position independent: 



Note that, for this to be possible, A must be negative, i.e., the potential is attractive. If one starts 
with a repulsive <5-potential, there is no way to deal with the divergent terms in the perturbative 
expansion of G. This reflects the fact that in dimensions D > 2 a repulsive ^-potential expels the 
S- waves from the Hilbert space Q. In this paper we consider only the attractive case. 
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n = -i hm km / — e e — > ; • (15) 

a ^ + £ ^ + J 2ix 2n ^ uj - (n + 1/2) eB + \i y ' 

Because of tke exponential in front of tke sum, one can close tke contour depicted in Fig. 1 
witk a semicircle of infinite radius in tke upper kalf-plane, and use residues to evaluate tke 
integral. Tke result, after taking tke limits e —>■ + and a — > + (in tkis order), is 

eB 00 

*o= -(n + 1/2) eB), (16) 

Z7V n=0 

wkere 8(x) is tke Heaviside step function. 

Witk respect to tke 'perturbed' part of n(x), it is easier to compute N± = J d 2 xni(x), 
wkere ni(x) = n(x) — n : 

Nl = -i km [ *± [ Sx G°(">*>°)Go("-M. (17) 

\ R 2-n ^ \2 eB J 

Performing tke integral over x witk tke kelp of tke identity J °° e~ z L m (z) L n (z) dz = 5 myn 
one finds: 

W/ (l — i£±i£\ 

N x = 1— km / due™ 12 , eB > N , (18) 

A<K 2 eB e-+o+J JL _ X ah l _ V v ; 

X R 2tt r \2 eB J 

wkere ip'(z) = dip(z)/dz. Tke integration over uj can also be performed using residues, but 
now tkere are two classes of poles to consider. Tke poles of tke first class kave tke form 
uj^ = —fi + (n + 1/2) eB (n = 0, 1, 2, . . .). Tkey are second order poles of ip', but are also 
simple poles of ip, and so are simple poles of tke integrand. Tkeir contribution to Ni reads 

oo 

N[ 1] = (n + l/2)eB). (19) 

n=0 

Tke poles of tke second class are given by tke roots of tke equation 

X R 2vr y V2 eB J K ' 

Examining tke grapk of tke digamma function |J one realizes tkat tke roots of Eq. ( |20| ) 
kave tke form uj® = — /x + (k n + 1/2) eB, wkere k < and n — 1 < k n < n (n — 1, 2, . . .). 
Tkeir contribution to Ni kas tke opposite sign: 

oo 

= Y t e( f i-(k n + l/2)eB). (21) 

n=0 

Taking into account an area factor A, one finally obtains tke following result for iV 
| d 2 xn(x): 

f eB \ oo oo 

N = ( — A - 1 ) % - {n + 1/2) eB) + £#(//- (k n + 1/2) eB). (22) 

^ ' n=0 n=0 



Tkis result kas a very simple pkysical interpretation: since tke potential kas zero-range, 
only tke S- waves are affected by it. Tkey are expelled from tke Landau levels (wkick kave 
eB/2n states per unit area), and mix among tkemselves to give new states witk energies 
equal to (k n + 1/2) eB. Tke explicit form of tkeir wave functions are obtained in tke next 
subsection. 
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B. Localized states 



The wave functions of the localized states can also be obtained from the Green's function. 
According to (§) and fllTf), 

^(x') — li m 2) — ^F^) G(lo; x, x') 

Go(a;f;x,0)G'o(^ 2) ;0,x / ). (23) 



(2) 



It follows from this and Eq. (|9|) that 



y 2-Kip'{-k e ) ^ Q k t -n 

For the lowest energy bound state (£ = 0), the generating function of Laguerre polynomials, 
Eq. (|A4]) , allows us to express the sum in ([24]) as an integral: since k < 0, we may write 



oo j I c \ oo roo 

„=o fc o - n ^ 7o 



/■OO ^ 

exp ( A; s — — - ) . (25) 



n=0 

00 ds 
Iq 1 - e" ; 

Changing the variable of integration from s to 2 = (e s — 1)~ , we find 

E = - / e-* z"^- 1 (1 + z) fco = -r(-fo) (-to, 1, € 

^ k -n Jo 

where f/(a, 6, z) is the Rummer's function which is singular at the origin || 



Combining (|J) and (p6|) we finally find 



^o(x) = -\J 2n ^_ ko) n-ko) M(x, 0) e" eBx2 /4 f/( _ fco; 1; ejBx 2 /2) . (27) 

This has precisely the form found by Perez and Coutinho [|J using the method discussed in 
the Introduction. 

Although the condition fco < was essential for obtaining (p6|) , this expression can be 
analitically continued for k (k 7^ 1,2, . . .), thus generalizing (P7| ) for the other localized 
states. Note also that, although they diverge at the origin (U(a, 1, z) ~ — [hiz +i/}(a)]/T(a) 
when \z\ — > 0), the wave functions ^(x) are normalizable. 
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C. Hall conductivity 



Now, let us consider the Hall conductivity. In the linear response approximation, it is 
given by |7j 

°2i{x) = e - <; hm o (D X2 -£>*,) J d 3 yG(x,y) yi G(y,x') \^ 

2 i 

= \ £ lim J J£ (D X2 - D* x ,) J d 2 y G(u; x, y) Vl G(u; y, x') |^ , (28) 

where D = V — ieA is the gauge covariant derivative, and D* is its complex conjugate. 
After performing the derivatives in x and the integration over y (for the latter, it is useful to 
use the identities [|[ L n (z) = L l n (z) — L^_ 1 (z) and ze~ z L x m {z) L\{z) dz = {n + 1) 5 m ,n), 
the 'unperturbed' piece of the Hall conductivity, obtained by replacing G with Go in (p8|) , 
reads 

4°i 0) = ^ Km lim /^e^5:[(2n+l)/ n 2 -2(n + l)/ n / n+1 ], (29) 



where 

e~ an 

fn = w- (n+ 1/2) eS + z/ (30) 

Performing the remaining integral (along the contour of Figure 1), and taking the limits 
e, a — > + , one finally obtains 

2 oo 

4°i 0) = -7r E e (A*-(" + l/2)eB). (31) 

27F n=0 

The 'perturbed' piece of the Hall conductivity, obtained by replacing G with G\ in 
Eq. (|28|), is easily shown to be zero (the integrand is an odd function of yi). One can also 
compute exactly the space average of the 'perturbed-unperturbed' pieces (in which one of 
the G's in Eq. (|2^) is replaced by Go and the other by G\)\ the calculation is rather tedious 
(it is sketched in the Appendix), but the result is remarkably simple: it is zero. 

Since all that remains is the 'unperturbed' piece of o"2i, we recover the remarkable result 
of Prange that the Hall conductivity of a two-dimensional electron gas is not affected by a 
5-impurity, even though such an impurity is capable of producing localized states. 

As a final remark, we notice that there is a simple reason why the correction to the Hall 
conductivity should vanish: as shown by Stfeda ||10|| , when the chemical potential is in an 



energy gap the Hall conductivity is given by following expression: 

ec 9N 

° u ~ ~~A ~dB' (32) 

where iV is the number of states below the chemical potential \l and A is the area of the 
system. With iV given by Eq. fl2"2"|) (and remembering that c = 1 in our units), it follows 
from the above expression and from Eq. ([H]) that a 2 \ = vfi^ ■ 
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Note added. The Green's function of a particle in a uniform magnetic field in the presence 
of a short-range impurity was previously obtained by Gesztezy et al (Gesztezy F, Holden 
H and Seba P 1989 On point interactions in magnetic field systems Schrddinger Operators, 
Standard and Non-standard ed P Exner and P Seba (Singapore: World Scientific) pp 147- 
164). We thank Prof. Pavel Exner for calling our attention to that paper. 



ACKNOWLEDGMENTS 



This work had financial support from CNPq, CAPES, FINEP and FUJB/UFRJ. 



APPENDIX A: 

In this appendix we sketch the calculation of the spatial average of cr^fx), obtained by 
replacing the first G in Eq. ( |28|) with Go and the second with G%. An explicit calculation 
shows that 



D x 2 ~ D l> 2 ) J d2 yG (uj;x,y) y x Gx(u);y,-x! 
eB yi G (y,0) 



Ajj 2vr r \2 eB 



x — x 



2 2 

+ {x 2 - Vz) £o(x, y) G? (0, x) - x 2 G (x, y) g (0, x) \, 



G (x,y) G (0,x) 



(Al) 



where Q is obtained from G by replacing the Laguerre polynomials L n (z) by their deriva- 
tives with respect to the argument, L' n (z). It follows from the explicit form of G and Q 
and from the identity L' n (z) = —L^^z) that, when calculating a?! = A 1 / d 2 xa^i (x), 



one has to deal with integrals which, after the change of variables (x, y) — > yeB/2 (x, y) is 
performed, have the following form: 

l£Z[P(*i,Vj)) = J d 2 x rfV- FM ^,%)i?(x 2 )i:((x-y) 2 )iI(y 2 ), (A2) 

where 



F(x, y) = i e^-Xj^ + - [x 2 + (x - y) 2 + y 2 ] 

and P(xi, yj) is a polynomial (of second degree) in Xi, yj. 

With the help of the generating function of Laguerre polynomials 



-l-Q 



exp 



tz 



Y,K{*)t n {\t\<i), 



n=0 



(A3) 



(A4) 



we define another generating function: 
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£ u m v n 



£,m,n=0 



xexp |_rf + (|t|iKW<1) . (A5) 



Performing the integrals over x and y, we obtain 



47i 2 rje v( - ap2+bci2+cp ' cl ~ iei: ' Piq ^ 
Z aPl (t,u,v; P , q ) = {1 _ ty+a{1 _ u y +P{1 _ v y^ (A6) 



where 



1 — uv 



(l-tt)(l-«) ; 

1 - ttt 

(l-*)(l-«)' 



(A7a) 
(A7b) 
(A7c) 



' 4(1 -tuv) v ; 

The integrals in ( |A2|) can then be obtained as the coefficients of the expansion in a power 
series in t, u and v of 

U i V ] P) QJ l p =q=0- 

(A8) 

The first type of integrals we need to evaluate is I^ n [— ixyyi + (iy\ + y\y2)/2\. Following 
the recipe given above, we obtain 

oo 



£,m,n=0 



d 2 i d 2 1 d 2 
dpidqi 

ATi 2 rf (-ic + ib) iir 2 (t - u)(l - v) 



' >^ + 2^ + 29A Zooo(t ' U ' Ui P ' q) "p— 



2 



(l-t)(l-u)(l-v) 4(1 -tuv) 

^— Y,{k + l)t k u k v k . (A9) 

4 fc=i 



It follows that 



in 2 

Q-^iyi + (ij/? + yiy 2 )/2] = — [{m + 1) fy m+1 (5 n , m - 5 n , m+1 ) - (m <-> £)}, (A10) 



so that (/ fc = [lu - (k + 1/2) eB + /j]- 1 ) 
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J2 lZ°nl-™ 1 y 1 + (iy 2 1 +ym)/2]f £ f m f n = 0. (All) 

£,m,n=Q 

The other two types of integrals we need can be found in an analogous way: 

I° mn [x 2 yi - ym] = — (m + l){Si >m - Si ;m+1 )(5 n>m - 5 n>m+1 ), (A12) 

T 



4 

Cnh^yi] = {£ + 1)(<W - <W+i)(<W - <W)- ( A13 ) 



Therefore, 

oo 

E (li'n-ij^yi - yM + lf\ m , n \-x 2 y^ f t f m f n = 0. (A14) 



£,m,n=0 



It follows from flXT]), (TATT|) , ( [AH] ) and the explicit form of G and G that 4i = 0. An 
analogous calculation shows that = 0, too. 
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Figure Captions: 

Figure 1: Integration contour in the complex w-plane used in the definition of the Feynman 
propagator. 
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